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Abstract 



The influence of charging effects on time-dependent transport in small 
semiconductor quantum dots with arbitrary level spectra is studied. Start- 
ing from an explicit time-dependent tunneling Hamiltonian, a non-Markovian 
Master equation is derived which is also valid in the nonlinear response regime. 
The many-body nonequilibrium distribution functions of the dot are calcu- 
lated and the I-V characteristic of the structure including the displacement 
currents is obtained. New resonant features show up in the Coulomb oscilla- 
tions and in the Coulomb staircase, and a new possibility to realize electronic 
pumps is described. 
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The influence of the Coulomb interaction on low-temperature quantum transport through 
metallic or semiconducting islands (quantum dots) has been the subject of many theoretical 
and experimental papers (see, e.g., and references therein). It shows up in a variety of 
interesting effects like, e.g., Coulomb blockade and resonant tunneling phenomena. However, 
the investigation of time-dependent perturbations and of fluctuations of the electrodynamic 
environement has started only recently, either neglecting charging effects |2|-[TT|| or consid- 
ering the metallic case for time translational invariant systems [T2|,[l3|]. An explicit time 



dependence of the applied voltages introduces a new energy scale hu in the problem, and a 
multitude of new effects is expected to occur, some of them relevant to device applications 
such as high-frequency oscillators. 

In this letter, we consider an interacting quantum dot with an arbitrary level spectrum 
subject to an explicit time- dependent classical field. It can describe a periodic modulation 
of the Fermi energy in the leads (i.e. time-dependent bias voltages) or time- dependent 
perturbations for the quantum states in the dot. We are especially interested in the effects 
of the Coulomb interaction in the limit of low tunneling rates but finite level spacing, so 
that resonant tunneling with thermally broadened line shapes will occur. In contrast to 
rfs. [0,0], we calculate the complete time- dependent (many-body) distribution function of 



the dot using a non-Markovian Master equation approach. As a consequence we obtain new 
features in the Coulomb oscillations which can not be seen in the time-independent case 
or for continuous level spectra like in metals. Furthermore, for certain asymmetric level 
structures in the nonlinear response regime, a new mechanism to realize an electronic pump 
is described. 

An important consistency check for our formalism is the fact that the sum of all currents 
into the system (including the displacement currents) is conserved ||. We have included 
the displacement currents in our formalism within the Coulomb-blockade model to make 
our calculations more realistic. However, since we are mostly interested in the case of low 
capacitances (where charging effects are important), we do not obtain a significant influence 
of the displacement currents on the I-V characteristics. 



As a model for an interacting quantum dot coupled to two reservoirs by tunnel junctions 
with capacitances Cl and Cr, we will use the time-dependent tunneling Hamiltonian H(t) = 
H (t) + V(t) + H T . Here, 

a=L,R k I 

describes non-interacting electrons in the reservoirs {a} and in the dot. A time-dependent 
shift of the Fermi energy of the electrons in the leads has been included in e^t) = e° a + 
A a (t). The tunneling part is given by Ht = J2 a J2ki{Tki a ia c i + c - c -)j where Tj^ denotes the 
tunneling matrix element, and the Coulomb interaction is approximated by the Coulomb- 



blockade model 14-16 



V(t) =E c (N D + n (t)) 2 (2) 

corresponding to the equivalent circuit (SET-transistor) schematically drawn in Fig. [l|. 
The particle number in the dot is denoted by = X^cjcj, Ec = e 2 /2C is the charging 
energy with C = Cl + Cr + C g and en (t) = CiYl^) + CrVr^) + C g V g (t) accounts for the 
time-dependent voltages of the left and right reservoirs eVLm(t) = ^l/r + A-L/Rif) as well 
as a time dependent gate voltage eV g (t) = \i g + A g (t) applied to the quantum dot by the 
capacitance C g . 

The tunneling currents from the left /right reservoir into the quantum dot are given by 

d 



m (3) 

where N a = J2k a \ a a ka is the particle number in reservoir a and p(t) denotes the density 
matrix of the system at time t. The expectation value < iV Q > can only change due 
to tunneling processes, therefore, the capacitive currents Ql/r (Ql/r is the charge on the 
left/right tunnel junction, see Fig. [D are not included in (|3]) and have to be added separately. 
The displacement currents are fixed by electrostatic considerations. The total current is 



-^(A L - A R ) + ~j^(A L - A g ) 

. Cr + C g rTun Cl r Tun t A \ 



and Ir is obtained by interchanging L and R. The capacitive current in the gate electrode 
is denoted by I g . A direct consequence is the conservation of all currents Ir + II + I g = 0, 
as emphasized in recent works by Biittiker et al ||. For C g = 0, we obtain at each t 
= —Iji(t), but the tunneling currents ([|) do not have such a property since the charge 
Q(t) on the island is a function of time [I7||. Generally, Eq. (|3]) gives the total current only 



if it is used to calculate the time-averaged current |L(|. If the capacitances are large, the 
first two terms on the r.h.s. of @ can be more important than the tunneling part. For low 
capacitances, on the other hand, the Coulomb energy is high and charging effects like the 
Coulomb-blockade will appear. We are interested in this limit and we will have to include 
the interaction (Q) in a nonperturbative way. 

In order to calculate the tunneling currents (||), we perform a time-dependent unitary 
transformation defined by H(t) = U{t)H(t)lJ- l {t) +ihU(t)U- x (t) with 

U(t) = exp(^ f dr[J2^a(r)N a + A D (r)N D ]) (5) 
ri J — oo a 

where A D (t) = C~ l Y,i=L,R, g CiAi{t) describes an effective time-dependent perturbation on 
the dot. The result is H(t) = H + V ' + Hrif) where Ho is the time-independent part of ([!]) (a 
constant shift of all energies has been included in ej), V = E c Nd(N d — 1) and H T {t) contains 
now time dependent tunneling matrix elements T^(t) = T^exp(| JL^ dr[A Q (r) — A^(r)]), 
and for an arbitrary interaction V on the dot the current lL/n{t) will depend only on the 
difference of the external voltages A^/^(t) — A ff (t). The differences that appear in the 
Hamiltonian, A a — A D , are linear functionals of A Q — A g , vie. A L — A D = [(Cr + C g )(A L — 
A g ) — Cr(Ar — A g )]/C, and analogously for Ar — A g . Eq. (|j) is unchanged by the unitary 
transformation since Nl/r commutes with U(t). A systematic perturbation expansion in 
Ht then yields 

C n W = eE f dtT rr ,(t,t')[N a (r') - N a (r)]P r (t') (6) 

rr 

where \r > are the eigenstates of H + V with energy E r , P r {t) =< r\p(t)\r > are the diagonal 
elements of the density matrix, N a (r) =< r\N a \r > is the particle number in reservoir a for 
the state |r >, and T rr >(t,t') are the time-dependent matrix elements 
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r rr /(t,t') = ^Re{< r\H T (t)\r' > 
n 

x < r'\H T (t')\r > ex P C-(E r - E r , + i^)(t - t'))} (7) 

for the transition r' (at time t') — > r (at time t). Non-diagonal elements of p(t) are of higher 
order in compared with the diagonal part P r (t) and do not contribute to the lowest 
order result This is equivalent to neglecting initial correlations between the dot and the 
reservoirs. Furthermore, we have included a finite lifetime r = ^ in ([?]) due to tunneling, 
where T = J2 a To and T a = 2irJ2k \Tki\ 2 5(E — €k a ) is the tunneling rate into reservoir a (for 
simplicity, we assume T a to be independent of the state I and the energy). This is equivalent 
to the well-known inclusion of collisional broadening effects in kinetic equations |18|,|19] . 

In the same way, we expand P r (t) in powers of and obtain the following non- 
Markovian master equation 

Pr(t) = E f dt'T rr/ (t,t')[P r ,(t') - P r (t)]. (8) 

Eq. (^) and (^|) are the result of our work in general form. They allow us to treat time- 
dependent, non-linear transport problems in the Coulomb blockade regime. Although they 
have been derived in the limit of low tunneling rates (r <C kT), the inclusion of a finite 
lifetime in (0) softens this restriction and leads to the exact result for the noninteracting 
case. In the presence of charging effects, it has recently been shown both experimentally 
22[ and theoretically J23| , that ^ describes the broadening of the Coulomb oscillation peaks 



in good approximation for discrete level spectra (AE > T) and temperatures (kT ~ T) that 
are not too low. For time- independent perturbations in the Markov- limit, we recover the 
well-known Pauli-Master-equation [f2(] . 

We will now approximate the interaction by the Coulomb-blockade model (Q), which has 
turned out to be sufficient to describe the charging effects in small quantum dots qualita- 
tively, provided that the number of electrons is not too small [|21|| . We therefore expect that 



it is also appropriate for the description of time- dependent phenomena. The eigenstates of 
/^o + V are then given trivially in the occupation number representation. 



For weak tunneling and large reservoirs, we can factorize P r {t) = Pl q P s (t) into an equi- 
librium part P? 9 for the reservoirs (described by the chemical potentials \x a according to the 
time-independent part Hq) and a part P s (t) for the dot. In order to find the stationary value 
for P s (t), we take a periodic modulation A a (t) = Aq sincut (A a = A a — Ad) and perform a 
Fourier transformation P s (t) = J2m=~oo Ps( m ) ex P( — imujt) with the result 

mhuj s ~ 



-)P s (m) = 5>( S| )JV 



m) 



+ 



L n al 



al 



) (2m(s) - l)(P s (n) +P S ») + (m -> -m) H 



(9) 



where s = |{^( s )}) denotes a state of the dot characterized by the occupation numbers rii of 
its single particle levels I, si is the state obtained from s by reversing the occupation of level 
I, i.e., ni(si) — 1 — ni(s), N s — J2i n>i(s) and E^j = ei + U^- — /i a + i^ where C/^ = 2NE C is the 
change of the Coulomb energy due to the addition of one particle to the dot. Furthermore, 
with (3 = 1/kT and f(E) = (exp f3E + l) -1 , we have defined 



P« m (P) = i m ~ n E J k+n (^-)J k+m (^-) Y(E + khu) 



Y(E) = -f(E) + — — 
v ; 2 JK J Am 



^ V 2 2vr ; 2tt ; 



(10) 
(11) 



where and tp denote the Bessel and digamma functions. Finally, the tunneling current 
reads 



fTun, 



m) 



-Y a V£N s P s {m) 



t-n l{s) ) Ps{n) + {m 



-m) 



n si 



(12) 



Eq. (^) is a linear set of equations for the Fourier components of the probability function 
P s {t). It can be solved numerically in a straightforward way and gives directly the time- 
dependence of all many-body distribution functions of the dot together with the tunneling 



currents (12). It can be applied to the nonlinear response regime (i.e. high bias voltages or 
high time-dependent perturbations), small or high frequencies, arbitrarily strong Coulomb- 
interaction and discrete level spectra. Furthermore, no assumptions are necessary for the 
distribution functions of the dot which can differ considerably from the equilibrium value. 
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For the time-independent case, Eq. @ reduces to the Pauli-Master equation which has been 
used in Refs. [p4|- p6| , pl|j to describe Coulomb oscillations and the Coulomb staircase in small 
quantum dots. The effect of time-dependent perturbations show up clearly in the argument 
E + k%uj (with E = Ep}) in flTop . It coincides with the one-particle excitation (tunneling) 
spectrum e; + + khu of the dot since the electrons can now absorb or emit energy quanta 
%uj. For the special case of one state in the dot, the solution of @ and fll2|) reads 



where A am = Fg m (E%) + F Q _ m It is identical to the result obtained in Ref. 



where interesting current oscillations have been derived from flllf ). Using our approach, 
we have studied this effect for many levels with Coulomb-interaction and have seen that it 
still remains observable for discrete level spectra (AE > bias voltage) or, for high charging 
energy, even in the case of continuous spectra like in metals. The capacitive currents will 
not change this picture significantly, since they are small for the low capacitances considered 
here. 

We will now discuss several applications of the formalism presented above. As a simple 
example we consider the DC-component of the current through a double-barrier quantum 
well with two degenerate states (ei = e 2 = e). Applying an AC gate voltage with frequency 
uj to the dot (A L — A R = 0, A g ^ 0), the Coulomb staircase reveals a fine-structure 
(Fig. 2) with additional steps at e + 2NE C + kTwo (k ^ ; N = 0, 1) due to the new 
excitations from the time-dependent perturbation. Even in the case of many levels in the 
dot with a continuous spectrum (where no fine structure due to the finite level spacing can 
be observed), these side steps will remain. In the same way, also the Coulomb oscillations 
will show additional resonances at a distance nhio (n = 0, ±1, ±2, . . .) from the conventional 
main peaks at e + 2NE C (N = 0, 1). However, the situation changes drastically if the two 
single particle levels are no longer degenerate (ei < e 2 ). In the static case, we get two 
Coulomb peaks at e\ and e 2 + 2E C (inset Fig. 2, dashed curve), i.e. the excitations e 2 and 
e\ + 2E C are hidden without time dependent perturbations. The reason is that on lowering 



the gate voltage (for small bias voltage and particle number N = 0), the first excitation 
will beiV = 0^iV = 1, occupying the level 1 = 1, i.e., the first peak is at e x . After this 
transition, the dot contains exactly one particle which can not escape due to the Coulomb 
blockade. The next possible excitation for iV = l^iV = 2isat energy 62 + 2E C and the 
ones with e<i and e% + 2E C are not visible. In the time-dependent case, the particle in the 
dot can absorb modulation quanta of energy hu and thus it can leave the dot or occupy 
the level / = 2. Therefore, the complete tunneling spectrum will show up in the Coulomb 
oscillations (see inset Fig. 2) in contrast to the static case. In other words, the effect of a 
finite charging energy in the presence of time dependent perturbations is not only a simple 
shift of the resonances at q (I = 1,2, . . .) by (I — 1)2E C but instead a splitting into many 
peaks at t\ + 2NE C (N = 0,1,2,.. .) where the amplitude is maximal at N = I — 1 and 
the peaks with N = I — 1 + k correspond to absorption (emission) of energy quanta khu. 
As in the case with the Coulomb staircase, these fine structures are not observable for a 
continuous spectrum in the dot where only the satellite peaks shifted by khu will remain. 

Another possibility is to apply an AC-signal to the reservoirs such that (At — Ap ^ 
0, A# — Ad = 0). Even for zero DC-bias voltage, we obtain a positive or negative DC- 
current through the quantum well due to the absorption or emission of energy quanta in 
the left reservoir. The effect is possible only if there is an asymmetry of level (s) in the dot 
with respect to the chemical potential (p = fiz = of the reservoirs: if the level closest 
to n is above /i, there will be a positive current, and vice versa. Consequently, the system 
acts like an electron pump in this case and the form of the Coulomb-oscillation peaks are 
changed dramatically as illustrated in Fig. |^. This example of a non-linear response is a 
generic result as can be seen by taking the (m = 0)-component of ([13]) for the case of a 
single nondegenerate level in the dot and expanding to second order in A a . For to — > 0, we 
obtain if c ~ /"( e Q — A i )(^ L ) 2 - The second derivative of the Fermi function will produce 
the positive and negative regions seen in Fig. |^. 

In conclusion, we have presented a formalism that allows us to calculate the frequency- 
dependent I — V^-characteristics of strongly interacting quantum dots. It takes into account 
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non-equilibrium many body distribution functions of the dot, it is not restricted to the linear 
response regime, and it is exact in the limit of small tunneling rates. We have also included 
the contribution of the displacement currents caused by the capacities and gate electrodes. 
The presence of five energy scales leads to a multitude of physical situations that will be 
studied in future work. 
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FIGURES 

FIG. 1. Equivalent circuit of the double-barrier structure in the static case. Cl,Cr are the 
capacitances between the dot and the leads, and C g is the capacitance between the dot and the 
gate electrode. Q is the excess charge on the dot, and the Qi are the charges at the capacitor 
plates. 

FIG. 2. DC current versus eF B ias for T = 5, E c = 75, w = 50, A L = A R = (dashed line), 
A L = A R = 50 (solid line), Cl = Cr = 10 -4 . The Coulomb staircase is modified by additional 
steps. Energies are measured in units of the total tunneling rate, r = Tl + Tr where Fl = Tr, 
currents in units of eT/h, frequencies in units of T/h, and capacities in units of e 2 /T. Inset: DC 
current versus eV g C/C g for two nondegenerate levels with ei ; 2 = ^12.5, other parameters as before. 
The Coulomb peaks are seen to develop side bands shifted by multiples of Tav. In addition, there 
are peaks corresponding to the single particle levels in the dot. 

FIG. 3. DC current versus eV g C/C g at zero DC bias voltage and A L = 30, A R = 0. Here, 
T = 2, uj = 20, and the other parameters are as before. The structure acts as an electron pump. 
The direction of the current depends on the position of the two many-particle levels with respect 
to the chemical potentials of the left and right leads. 
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